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Abstract: We analyze the Type-II two Higgs doublets model in light of the newly 
discovered Higgs-like particle with mass 125 GeV. The observed 125 GeV particle is 
^ identified with the light CP-even Higgs boson in the two Higgs doublets model. We 

Q study the parameter space of the model consistent with the Higgs data, branching 

■^ ratio of -B — 7- Xs7 as well as precision electroweak measurements. We also incorporate 

^ theoretical constraints — perturbativity of the couplings and vacuum stability, in our 

^ — I study. We find that only a small parameter space of the model remains viable. The 

^ phenomenology of the heavy Higgs bosons in the surviving parameter space is studied. 
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1 Introduction 

Recently the ATLAS and CMS collaborations announced the discovery of a Higgs-like 
particle with a mass M^ — 125 GeV [1, 2]. Evidence of this new particle has also 
been reported by the CDF and D0 collaborations [3]. However, it is far from certain 
that this newly discovered particle is the standard model (SM) Higgs boson responsible 
for electroweak symmetry breaking. The couplings of this Higgs-like particle deviate 
(although not statistically significantly) from SM expectations [4, 5]. Even with the 
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updated measurements [6-13], the deviations still remain. Thus it is possible that this 
Higgs-like particle is a hint of new physics beyond the standard model. 

A particularly well motivated class of new physics is the two Higgs doublet model 
in which electroweak symmetry is broken by two elementary scalar fields. Famously, 
the minimal supersymmetric standard model (MSSM) contains two higgs doublets to 
account for masses of all quarks and leptons, and the parameters are constrained by 
supersymmetry. However, in this work we consider a more generic two Higgs doublet 
model. For a recent review of a general two Higgs doublet model see Ref. [14] and 
references therein. In particular we will focus on the CP conserving type II two Higgs 
doublet model (2HDM-II) in which one scalar field couples only to the up-type quarks 
and the other couples to the down-type quarks and leptons. 

There exists a large literature on the properties of the neutral scalar boson couplings 
in light of the Higgs data; see for example Ref. [15-18]. In this work, in addition to 
performing a global fit to the currently available Higgs data, we also consider the viable 
parameter space of the model. Specifically, we assume that there is no other states 
except those of the 2HDM-II up to some cutoff scale, A. Thus if there is no viable 
parameter space for a specific value of the cutoff, we can conclude that, if there is 
new physics beyond the standard model, the 2HDM-II cannot be the only new physics 
below that cutoff scale. 

From theoretical view point, the model has to allow for an electroweak symmetry 
breaking vacuum. We also impose a constraint on perturbativity of the coupling con- 
stants of the model. What we mean by perturbativity will be made clear in section 4.1. 
We insist that perturbativity must be satisfied at all energy scales up to the cutoff 
scale [19]. This is different from Ref [17] which seems to impose perturbativity only at 
the electroweak scale. 

Existing experimental data also constraint possible new physics. The absence of 
large flavor changing neutral interactions places a strong bound on the mass of the 
charged Higgs boson. Similarly, the success of the standard model in describing preci- 
sion electroweak measurements constraints possible new physics states. In this work, 
we will utilize both experimental and theoretical constraints in determining a viable 
parameter space for the 2HDM-II. 

The paper is organized as follows. In section 2 we briefly describe the 2HDM- 
II and we set our conventions and notation. Then we perform a global fit to the 
Higgs data in section 3. In section 4 we subject the 2HDM-II to both theoretical and 
experimental constraints to determine the viable parameter space. We briefly discuss 
the phenomenology of the heavy CP-even neutral Higgs boson, H, in the surviving 
parameter space in section 5. We then conclude in section 6. 
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2 The Model 



Here we give a brief overview of the 2HDM-II and we set our notation. We take the 
two scalar doublets $i and $2 to have hypercharge 1/2. They can be expanded as 






+ Pj + ir]j)/V2^ 



The Yukawa coupling to fermions are given by 

^yuk = -ynQhiicr'^^DuR - VdqL^ldR - VeGL^ieR + h.C. 

Expanding out the CP-even neutral scalar sector we obtain 
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where tan/3 = ^2/^1 and f^ = f^ + f|. The two CP-even neutral scalars mix with each 
other. The mass eigenstates are given by 

cosa sinaA /H^ 
— sin a cos a J \ h 
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where h is the lighter eigenstate to be identified with the observed 125 GeV Higgs-like 
particle. In terms of these mass eigenstates, we find 
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The couplings of these two eigenstate h and H to gauge bosons can be obtained 



from the kinetic terms of $1 and $2- They are 
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2.1 Scalar Sector 



The scalar sector of the 2HDM is the most model dependent part. Here we will focus 
on the simplest scalar potential consistent with CP symmetry 
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It is more convenient to characterize the scalar sector by their physical masses and 
the mixing angles — Mh, Mh, Ma, Mh±, a, tan/3: 

■> 

sin a Ml, + cos aMrr = ^-7:^1, 

" 1 + tan^ /3 



tan^ B 
'^^"^^""^^=l + tan2/3""^' 



cos^ a Mf + sin^ a Mi = -tt^v^X 



{Ml — Mjj) cos a sin a = (A3 + A4 + A5)- y^v 

Ml = -\,v\ 



tan/3 2 (2- 

1 + tan^ /3 



M^± = -^(A4 + \,y = Ml + ^(As - \,)v\ 

The set of parameters we use in our parameter-space scan consists of a, tan/3, M^, 
Mh, Ma and Mh±. We identify the light CP-even scalar, /i, with the observed 125 
GeV resonance. Hence, for the remaining part of this work we set Mh = 125 GeV. 

3 Fit to the Higgs Data 

The mass parameters M^, Ma and Mu± affect Higgs data observables only through 
suppressed radiative corrections. We determine the VEV ratio tan/3 and the neutral 
scalar mixing angle a from all the reported Higgs data. 

Experimental data are reported in terms of a signal strength, /i, defined as 

— ^ ^^ 

When the signal strength is not directly reported by the experimental collaboration, 
we extract it from the reported 95% exclusion limit following the procedure given in 
Ref . [20] . Here we briefly review the procedure and refer the reader to the reference for 
details. The experiments report the upperbound on the rate at 95% C.L., -Robs, and 
the expected upperbound at 95% C.L. in the absence of the Higgs boson, Rexp- The 
signal strength and its uncertainty, a, are given by 

/^ = Robs — Rexp, 0" = -r~F:7^- (3-2) 

1.9d 

We collect the signal strengths for each search channel in Table 1. 

3.1 Production Cross-sections and Branching Ratios 

Here we work out the Higgs production cross-section and branching ratios for a non- 
stadard model Higgs coupling. At tree-level, the non-standard Higgs couplings to SM 
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Channel 


Signal Strength (//) 


Reference 


ATLAS 77, 7 TeV 


1 6+°-^ 

-■-•0-0.8 


[21] 


CMS 77, dijet-tagged, 7 TeV 


4.21 ±2.04 


[22] 


CMS 77, untagged 0, 7 TeV 


3.15 ±L82 


[22] 


CMS 77, untagged 1, 7 TeV 


0.66 ±0.95 


[22] 


CMS 77, untagged 2, 7 TeV 


0.73 ±1.15 


[22] 


CMS 77, untagged 3, 7 TeV 


1.53 ±1.61 


[22] 


CMS 77, dijet-tight, 8 TeV 


1.32 ±1.57 


[22] 


CMS 77, dijet-loose, 8 TeV 


-0.61 ±2.03 


[22] 


CMS 77, untagged 0, 8 TeV 


1.46 ±1.24 


[22] 


CMS 77, untagged 1, 8 TeV 


1.51 ±1.03 


[22] 


CMS 77, untagged 2, 8 TeV 


0.95 ±1.15 


[22] 


CMS 77, untagged 3, 8 TeV 


3.78 ±1.77 


[22] 


ATLAS 77, 8 TeV 


1.6 ±0.32 


[11] 


ATLAS ZZ, 7 TeV 


1 4+1.3 
-■-■^-0.8 


[23] 


CMS ZZ, 7 TeV 


6+°-^ 


[24] 


ATLAS ZZ, combine 7 & 8 TeV 


1.5 ±0.6 


[11] 


CMS ZZ, combine 7 & 8 TeV 


qi+0.30 

U.yi_Q24 


[12] 


ATLAS WW, 7 TeV 


n 1 +o-'^ 

^■-■—0.6 


[25] 


CMS WW, 7 TeV 


0.4 ±0.6 


[24] 


ATLAS WW, 8 TeV 


1.45 ±0.56 


[8] 


CMS WW, combine 7 & 8 TeV 


0.76 ±0.21 


[13] 


CMS bb, 7 TeV 


1 9+2-1 
-■-•^-1.7 


[24] 


CMS bb, 8 TeV 


1.07 ±0.62 


[9] 


Tevatron bb 


2.0 ±0.7 


[10] 


ATLAS Tf, 8 TeV 


0.7 ±0.7 


[26] 


CMS Tf, 8 TeV 


0.88 ±0.50 


[9] 



Table 1: The signal strengths and the corresponding error for the Higgs data used in 
the fit. 



fields can be characterized by rescaling coefficients, q's, relative to the standard model 
higgs couplings as follows 

h , 2 -\- ^ 2 \ h - h - h _ h 

Ch = cv- ylmy^W^W^^ ± nizZnZ^) - Ct-nittt - Cb-nibbb - Cc-rricCC - Cr-rrirTT 

V /^^ ' V V V V 



(3.3) 
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In this work we assume that other scalars are sufficiently heavy that their effects on 
the 125 GeV Higgs boson decay channels are negligible. 

The main production channels considered here are gluon fusion (ggF), vector boson 
fusion (VBF), vector boson associated production (Vh) and tih production. These 
production cross-sections are given in terms of the SM ones by 
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The rescaling factors for the partial decays widths are 



bb 



\^b\ ^ bb ' 



99 






T^SM 

gg 



1.03Q - 0.05C6 



cc ' 
2 



|„ \2-pSM p 



vv 



l^vl J- vv^ 



\Ct\ 



^SM 



(3.4) 



1.03-0.05 



±ji_ 
-rSM 

77 



1.03ci - IMcv 



§1.03 



1.04 



(3.5) 



where we have assumed that the loop-induced, charged Higgs contribution to V^^ is 
negligible. This assumption can be justified once we include experimental constraints 
from section 4.2. The rescaling factors for the 2HDM-II are given by 



cos a 



sma 



Q 



Cb 



I ) 



cy = sin(a + (3). 



(3.6) 



sin /3 cos (3 

For completeness we include the rescaling of the partial decay width for h — )■ 7Z 



rSM 



-(0.38 + 0.37z)0.86q - (3.32 + 3.98i)cv 



(3.7) 



-(0.38 + 0.37i)0.86 - (3.32 + 3.98z) 

The general expressions for the loop-induced decay of a neutral scalar boson are col- 
lected in appendix A. 

3.2 Global Fit 

We determine the best ffi values of tan /3 and a using all the reported Higgs data from 
the Tevatron, LHC 7 TeV and LHC 8 TeV runs. The theoretical SM Higgs boson 
predictions for the cross-sections and branching ratios are taken from the LHC Higgs 
Cross-section Working Group [27]. We obtain best ffi values tan/3 ~ 0.01 and a ~ 1.56. 
At this small value of tan /3 the top Yukawa coupling is non-perturbative at the weak 
scale. Thus the best fit value doesn't seem to correspond to viable model parameters. 
Assuming a Gaussian distribution around the best fit values, we determine the 
region of parameter-space consistent within 1- and 2-a with the best fit values. The 
result is shown in figure 1. The dashed line indicates the decoupling limit, a + P = 
71 /2, where the couplings are SM-like (up to a possible sign fiip for the down Yukawa 
couplings). The figure is truncated at a < 0.5, since for larger values tan/3 decreases, 
forcing the top Yukawa coupling into a non-perturbative regime. 
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^ 




Figure 1: The region of parameter space witliin 1- and 2-0" of tlie best fit values. Tlie 
dashed line is the decoupling limit, a + (3 = ±7r/2, where the couplings are SM-like (up 
to a possible sign flip for the down Yukawa couplings). 



4 Bounds on Parameter Space 

4.1 Theoretical Constraints 

We impose constraints amounting to the potential being bounded from below. To this 
end the couphngs must satisfy [19] 



Ai,A2 > 0, 

A3 > -V A1A2 , 

A3 + A4 - IA5I > - V A1A2 , 

at all scales up to the cutoff scale A. 

We also impose the following perturbativity constraint on the couplings 



(4.1) 



47r 



< 1 



4:71 



< 1. 



(4.2) 
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We insist on these constraints up to the cutoff scale for all the Yukawa and scalar 
couplings. We list the beta-functions used in evolving the coupling constants in ap- 
pendix B. 

4.2 Experimental Bounds 

A wealth of experimental data, particularly from precision measurements, places strong 
constraints on the spectrum of the 2HDM-II. A newly published result on a direct search 
for the charged Higgs at LEP yields the 95% CL lower bound Mh± > 80 GeV [28]. At 
present there is no lower bound on the charged Higgs mass from the Tevatron or LHC 
data. A much tighter constraint on the charged Higgs mass can be deduced from rare 
decay processes. By analyzing the branching ratio Br{B — )■ ^^7), Ref. [29] obtained 
the bound Mh± > 380 GeV at 95% confidence level. A direct search at LEP also places 
a 95% limit Ma > 93 GeV [30]. 

Electroweak precision measurements also place strong constraints on the spectrum 
of a 2HDM. We concentrate on the oblique S and T parameters. In the Standard 
Model, the best fit values for the S and T parameters for M^ ~ 125 GeV and Mf = 173 
GeV, as well as their correlation matrix {Mcorr) are [31]. 

5 = 0.03 ±0.10, ^ ^/l 0.891\ 

T = 0.05 ±0.12, "" VO-891 1 / 

Additional contribution to the S and T parameters from the heavy Higgs bosons are 
given in Ref. [14]. 

4.3 Viable Parameter Space 

In this section we present the result of our parameter space scan consistent with both 
the Higgs data and the experimental and theoretical constraints discussed above. A 
point in parameter space is a set of values of the parameters Mh, Ma, Mh±, a, tan/3, 
as discussed in section 2.1. For every parameter point consistent with the experimental 
constraints of 4.2, we determine the corresponding scalar and Yukawa couplings using 
equations (2.5) and (2.8). These couplings are then evolved numerically from the weak 
scale, u^, to the cut-off scale. A, using beta-functions listed in appendix B. Finally the 
couplings are checked against the theoretical constraints of section 4.1. 

The interplay between the experimental and theoretical bounds on the spectrum 
can be easily understood. Intuitively, experimental bounds — Br{B — )■ Xg'y) and elec- 
troweak precision data tend to drive the mass of the Higgs bosons heavy in order to 
leave a small imprint on low energy observables. Moreover, a non observation of neu- 
tral scalars at LEP means Ma ^ 93 GeV. At the same time, the more massive the 
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Figure 2: The viable particle spectrum for 0.75 < tan/3 < 2.75. The spectrum depends 
non-trivially on the charged Higgs mass (Mj^i). In each plot the viable parameter space 
in Mh-Ma plane shrinks as Mh± increases. Mh± = 380 (420) GeV corresponds to the 
color blue (red). 



spectrum is, the larger the scalar couplings. Perturbativity constraints limit how large 
these scalar couplings can be, hence limit from above the Higgs spectrum of the theory. 
Combining these experimental and theoretical considerations, we find that with the 
cut-off scale (A) at 2 TeV, there is no viable parameter space. Reducing A to 1 TeV 
opens up a small parameter space for small values of tan/3. Thus, for the rest of this 
work we will take A = 1 TeV. The viable spectrum is shown in figure 2. 

We end this section with a brief discussion of the sensitivity of the viable parameter 
space on the choice of perturbativity condition, eg. equation (4.2). Had we imposed 
instead that all the reduced couplings remain less than 1/2, there would be no viable 
parameter space. We find that relaxing the perturbativity constraint from 1/2 to 3/4 
opens up a small viable parameter space for the cases tan/? = 1, 1.25. Relaxing this 
constraint further to 1 leads us to the viable parameter space that we have in figure 2. 
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Figure 3: The heavy CP-even Higgs boson production cross-sections. For each plot, 
a is taken to be 0.78 (0.58) for tan/3 = 1.0 (1.5). This choice of a minimizes x^ltan^ 
for a given value of tan/3. 



5 Phenomenology Of The Other Higgs Bosons 

In this section we study the phenomenology of the neutral CP-even and CP-odd Higgs 
bosons, H and A. We will focus on their production cross-sections and branching ratios. 
For a wide range of production and decay channels, we can deduce the cross-sections 
and branching ratios by scaling from the corresponding quantities for the SM Higgs 
boson. When this scaling procedure is not available, we compute the corresponding 
quantity at leading order. 

5.1 Phenomenology of the CP-even H 

The couplings of the heavy CP-even Higgs boson H to the SM fermions and gauge 
bosons are modified with respect to the SM Higgs boson couplings by the rescaling 
factors 
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Figure 4: The heavy CP-even Higgs boson branching ratios for the case of heavy Ma- 
For each plot, a is taken to be 0.78 (0.58) for tan/3 = 1.0 (1.5). This choice of a 
minimizes x^ltan/3 for a given value of tan/3. 
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Figure 5: The heavy CP-even Higgs boson branching ratios for the case of light Ma- 
Here a is taken to be 0.78. This choice of a minimizes x^|tan/3 for a given value of tan /3. 



The H production cross-sections can be readily obtained by rescaling the corresponding 
calculations for SM Higgs boson production. The CP-even H boson production cross- 
sections for the 8 and 14 TeV LHC are shown in figure 3. Note that since the mixing 
angles a and /3 are close to the decoupling limit, the H production cross-sections in the 
vector boson fusion and the W- and Z-associated production channels are suppressed. 
In a large portion of the available parameter space Mh is lighter than 2Mh, Ma 
and Mh±; see figure 2. When this is the case all the available decay channels for the 
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heavy Higgs boson, H, are the same as for the SM Higgs boson. The branching ratios 
can be obtained by a simple rescahng from the corresponding values of the SM Higgs 
boson. When Mh is sufficiently heavy, the decays H — )■ hh, H — )■ AZ and H — )■ AA 
become available. The partial widths are^ 



T{H -^ 



87271 Mh V Mf, 



T{H ^ AZ) = -^ sm\a + /3) 

8v27r 



2 

2^ 
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Ml 



3/2 



(5.2) 



where Gp is the Fermi constant, cj) = h, A and 

(2M2 + Ml) cos(a + (3) sin(2a) 



^Hhh 
^HAA 



2 sin (3 COS (3 
{3Ml + 2Ml) sin(a - /3) + {Ml - 2Ml) sin(a + 3/3) 



(5.3) 



4 sin /3 cos /3 

The branching fractions for H decays are shown in figure 4 and 5. Since the mixing 
angles a and (3 are close to the decoupling limit, the decays into a pair of massive vector 
bosons are suppressed. When H ^ A + X is kinematically forbidden, the H boson 
decays predominantly into quarks and gluons; see figure 4. In this case, it is most 
likely to observe the i7-boson in di-photon decays. However, when H -^ A + X is 
kinematically allowed, the decay mode H — )■ AZ becomes the most dominant; see 
figure 5. As we will show in the next section, the A has sizable branching ratios into 
'-yZ and 77. Thus it might be possible to observe H — )■ AZ in a photon plus four 
leptons or two photons and two leptons channels. 

5.2 Phenomenology of the CP-odd Boson 

The pseudoscalar A does not couple at tree level to a pair of electroweak gauge bosons. 
Thus its main production mechanism is via gluon fusion. Its production cross-section 
can be obtained from that of the Higgs in the Standard Model by rescaling by a factor 

2 



' g — 



cot/?^"^'"> 



where r = M\/4M^ and the function /(r) is defined in appendix A. The production 
cross-sections of the pseudoscalar A at the LHC with 8 TeV and 14 TeV center-of-mass 
energy are shown in figure 6. 



^Our result does not agree with that of Ref [32] whose partial width for ■!> — )■ (j)Z does not have a 
correct mass dimension. 



-12- 



tan/?= 1.0 




tan/3= 1.5 



100 150 200 250 300 350 400 450 
MA[GeV] 




300 350 

M^[GeV] 



Figure 6: The heavy CP-odd Higgs boson production cross-sections. For each plot, a 



is taken to be 0.78 (0.58) for tan/3 = 1.0 (1.5). This choice of a minimizes x^ 
a given value of tan /3. 
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The CP-odd boson A has a mass in the range 93 < Ma ^ 470 GeV, see figure 2. A 
light A can decay into a pair of light quark-antiquark as well as into two photons, much 
like the SM Higgs. For sufficiently large Ma-, A can also decay into pair consisting of 
a neutral CP-even scalar boson and a vector boson, 0Z, where = /i, if . The partial 
decay widths for A — )■ gg is 



Y{A ^ qq) 



GpMj, , 



4^2 



TT 



AMVM] 



Ai 



(5.5) 



where 5q = cot/3 (tan/3) for up-type (down-type) quark. The partial decay width into 
(pZ is given by 

r(A -> 0Z) = -^52 



8v^7r 



{Ml - {Mz - M^f) {Ml - {Mz + M^)- 
Ml 



3/2 



(5.6) 



where (5/i {5h) = cos(a-|-/3) (sin(a-|-/3)). The loop-induced branching ratios into gg, 77 
and 7Z are given in appendix A. The CP-odd A branching ratios are shown in figure 7. 
Note that the suppression of Br {A — )■ hZ) arises because the two mixing angles a and 
/3 are close to the decoupling limit, see equation (5.6). The dominant decay channels 
for A are hh and gg for a low mass A and HZ and ti for a heavy A. 

It is also interesting to note that for the low mass range below the ti and HZ 
thresholds, the A production cross-section and branching ratios into 77 and 7Z are 
enhanced compared to the corresponding Standard Model Higgs boson counterparts. 
Our estimate of the signal strength in these two modes is shown in figure 8. With the 
large signal strength for Ma > 160 GeV, this scenario could be excluded using current 
Higgs data. However, at the moment both the ATLAS and CMS collaborations only 
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Figure 7: The CP-odd Higgs boson branching ratios. For each plot, a is taken to be 
0.78 (0.58) for tan/3 = 1.0 (1.5). This choice of a minimizes x^|tan,g for a given value of 
tan/3. Diagrams on the top correspond to the case when Mu is sufficiently heavy and 
cannot be in the decay product of A. Diagrams on the bottom are when Mu is light 
enough to be in the decay product. For definiteness, we take Mu = 130 GeV. Here we 
ignore the effect of near threshold production 



provide the 95% exclusion limit on the signal strength in these channels up to a mass 
of 150 GeV [33-36]. 

5.3 Phenomenology of the Charged Boson H"^ 

The viable mass of the charged Higgs boson is larger than that of the top-quark. Thus 
its main production cross-section is from gb — )■ H^i and gg — t- H^ib [37]. We use 
Madgraph to compute the lowest order production cross-sections for both processes. 
We follow Ref. [38] in keeping the renormalization scale and the PDF factorization scale 
fixed at fi = Mt + Mh± . For comparison, we also compute the production cross-section 
from qq — )■ H^ib where q is the valence quarks. The production cross-sections are 
shown in figure 9. 
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The physical charged Higgs couplings to quarks are given by 

{V2GFy/^ [cot /3{m^),UiV*jPLdj + tan /3{md)jU,V*j Pad j] H+ + h.c, (5.7) 

where m„ (fnii) is the mass matrix of the up (down) type quarks and V is the CKM 
matrix. Thus its partial width r(77+ — )■ tb) at leading order is 
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The charged Higgs can also decay into a neutral scalar and a vector. In this case the 
partial width can be read off from equation (5.6) with the replacement Ma — )■ Mh± 
and 5a = 1- The branching ratios of the charged Higgs boson are shown in figure 10. 
Much like for A — )■ hZ, the branching ratio for H^ — )■ hW^ is suppressed because a 
and (3 are close to the decoupling limit. For a light H{A)^ the decay H^ — )■ H{A)W'^ 
is dominant. This decay leads to an event with 3 or 4 6-quarks (depending on the 
production channel) plus a charged lepton and missing energy which could be searched 
for at the LHC. 
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Figure 9: The charged Higgs boson production cross-sections for the case tan/3 = 
1.0 (1.5) for the left (right) plot with a = 0.78 (0.58). a is chosen such that it 
minimizes x^ltan/3 for a given value of tan/3. Diagrams on the top correspond to the 8 
TeV LHC while the bottom are for 14 TeV. 



6 Conclusions and Discussions 

We have performed a parameter scan for the CP conserving 2HDM-II consistent with all 
the available Higgs data. We take into account theoretical bounds — vacuum stability 
and perturbativity of the couplings, as well as experimental bounds from electroweak 
precision measurements and Br{B — > X^'-)) in our scan. We use a working assumption 
that the 2HDM-H is a valid low energy effective theory up to a cut-off scale A. We 
found that for A > 2 TeV, there is no viable parameter space consistent with all the 
mentioned constraints. However, if we assume A = 1 TeV, only a small parameter 
space of the 2HDM-H is consistent with all the mentioned constraints. In particular, 
the ratio of the vacuum expectation value of the two scalar doublets, tan/^, lies in the 
range 0.75 < tan /3 < 2.75. 

Our results show that the charged Higgs boson is the most constrained sector of the 
2HDM-II. The perturbativity constraint demands Mii± < 420 GeV while constraint 
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Figure 10: The charged Higgs boson branching ratios for the case tan/3 = 1.0 (1.5) 
for the left (right) plot with a = 0.78 (0.58). a is chosen such that it minimizes x^|tan/3 
for a given value of tan/3. The masses of the neutral Higgs boson are chosen so that 
they are consistent with the viable parameter space shown in figure 2 



from Br{B — )■ ^^7) pushes Mh± > 380 GeV. Since the mass of the charged Higgs 
boson is much heavier than that of the top-quark, its main production mechanism is 
production in association with the top or in association in the top and the bottom- 
quarks. The charged Higgs, once produced, decays mostly into a top- and bottom- 
quark pair which makes its detection difficult. However, in a corner of parameter 
space where the heavy CP-even H (or the CP-odd A) is light enough, the decay into 
H{A)W becomes comparable to the tb channel. The light H{A) then mostly decays 
into a bb pair. This decay chain leads to an event with 3 6-quarks (or 4, depending on 
the production mechanism) plus a charged lepton and missing energy which makes it 
possible to be searched for at the LHC 

The neutral scalars sector is not as tightly constrained. The masses of both the 
heavy CP-even Higgs boson, H, and the CP-odd, A, can take on values in a large range; 
see figure 2. The production cross-sections for these two particle are sizable at both the 
8 TeV and the 14 TeV LHC. Since the mixing angles a and /3 are close to the decoupling 



-17- 



limit, the H is mainly produced via gluon fusion. Hence a low mass if, where H decays 
predominantly into hh, will be difficult to observe at the LHC. However, for a heavy H 
the decay H — )■ AZ becomes available and could lead to spectacular decay signatures 
of one photon and 4 leptons or two photon and leptons if the A decays into 7Z or 77 
and the Z's decay leptonically. 

A low mass pseudo scalar A has an enhanced production cross-section compared 
to the SM Higgs boson counterpart. Moreover, its branching ratios into 7Z and 77 are 
also enhanced due to the absence of massive gauge boson decay channels. This leads 
to a large signal strength for the 7-Z' and 77 channels. Therefore, current LHC data 
can exclude (or establish) this scenario. 

Lastly we emphasize that our framework assumes 2HDM-II to be an effective low 
energy theory of a more complete theory. Our analysis suggests the 2HDM-II descrip- 
tion of the electroweak symmetry breaking is only valid up to a scale around 1 TeV. 
Hence, if nature chooses to break electroweak symmetry by the 2HDM-II, there must 
be new particles waiting to be discovered with masses around 1 TeV. 
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A Loop-induced Decays of Neutral Scalar Bosons 

The expressions for loop-induced decays of the Standard Model Higgs boson, as well 
as for the MSSM neutral bosons decays, are given in Ref. [32, 42]. Here we give the 
rescaling factors for the partial decay width suitable for the 2HDM. 
The rescaling factors for — )■ gg^ 77 where (p = h, H are 

E/cf ^1/2(1-/) 



-^ 99 

rSM 
99 



-^ 77 
rSM 

77 



EfNcQ^fAi/^iTf) - 4A^iTw) + X^h+h-Ao{th±) 



EfNcQ}Ai/2{rf)-A,{Tw) 



(A.l) 



where Nc and Qf are the number of colors and the electric charge of fermion /. X^h+h- 
is the trilinear coupling of to the charged Higgs boson in units of 2M?/f . The 
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functions Ai give the 1-loop contribution of a spin-i particle. They are given by 
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with Ti = Ml/AMf ioii = f,W and H^. 

The rescahng factor for — )■ 7Z where (p = h, H is 






E/ iVc ^^^'^y ^^^^ cf A?^,(r^, A^) + 4At(rH^, A^) + ^A<^^.^-Ao^(r^±, A^^) 



E/ Ar^ Q/(^^.HQ/-^) At^^(^^, A,) + Atirw, \w) 



(A.4) 

where sw{cw) is the sine (cosine) of the weak mixing angle, /? is the left-hand weak 
isospin of the fermion /, r, = 4Mj^/M| and Aj = 4Mf /M|. The function Aj's are 
defined as 
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where the function Ji and I2 are defined as 
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The function / is the same as in the — )■ 77 case. The function g is defined as 
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For the CP-odd boson, its loop induced partial decay width into gg, 77 and 7Z 
are given by [32] 
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B Beta- functions 

Here we list the one-loop beta-function for the gauge, Yukawa and scalar couplings. 
For the Yukawa couplings, we consider only the third generation contributions. We 



use the shorthand notation at 
couphng beta-functions are 
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Finally the scalar couplings beta-functions are 
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